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Abstract 

We present an analytical description of the phase transitions from a nucleon gas to 
nuclear matter and from nuclear matter to quark matter within the same model. The 
equation of state for quark and nuclear matter is encoded in the effective potential of a 
linear sigma model. We exploit an exact differential equation for its dependence upon the 
chemical potential associated to conserved baryon number. An approximate solution 
for vanishing temperature is used to discuss possible phase transitions as the baryon 
density increases. For a nucleon gas and nuclear matter we find a substantial density 
enhancement as compared to quark models which neglect the confinement to baryons. 
The results point out that the latter models are not suitable to discuss the phase diagram 
at low temperature. 
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1 Introduction 



The equation of state for strongly interacting matter at nonzero density is needed for the 
understanding of neutron stars [1], as well as for the interpretation of heavy ion collision 
experiments [2]. Any analytical description of the equation of state for baryons at nonzero 
baryon density has to cope with the problem that the effective degrees of freedom change from 
nucleons at low density to quarks at high density. We attempt in this note a unified description 
of both the nuclear gas-liquid transition and the transition to quark matter. For this purpose 
we work within an effective linear meson model coupled to quarks and nucleons. It should 
describe the low momentum degrees of freedom of QCD for the range of temperatures and 
densities which are relevant for these phase transitions. Our main computational tool will be 
a new exact functional differential equation for the dependence of the effective action on the 
baryon chemical potential and an approximate solution to it. We will see many similarities but 
also important differences as compared to a mean field theory treatment. 

Our main interest are the chiral aspects of the equation of state for quark and nuclear 
matter at nonzero baryon density and the order of the involved phase transitions. In this note 
we concentrate on the simplest possibility where only the color-singlet chiral condensate plays 
a role in the transitions. Our investigation should constitute a useful point of comparison for 
other models with more complicated condensates. In fact, in addition to the nuclear and quark 
matter phases a number of interesting possibilities like the formation of meson condensates or 
strange quark matter [3, 4, 5] have been proposed. Also an extensive discussion has focused 
around the symmetry of the high density state, where the spontaneous breaking of color is 
associated with the phenomenon of color superconductivity [6, 7, 8, 9]. Other ideas concern 
the spontaneous breaking of the color symmetry in the vacuum [10, 11]. In this note we 
adopt the working assumption that possible additional condensates have only little influence 
on the transitions associated with the order parameter of chiral symmetry breaking 1 . We also 
concentrate mainly on the case of two quark flavors and neglect isospin violation 2 . 

Within this setting and using a rather crude approximation for the transition from quark 
to nucleon degrees of freedom, we find that for low temperature both the nuclear gas-liquid 
and the hadron-quark transitions are of first order, in accordance with indications from ear- 
lier investigations (cf. [13, 12, 14] and [9] and references therein). The phase transition from 
nuclear to quark matter tends to be much stronger (larger surface tension) than the gas-liquid 
nuclear transition. The first order character of these phase transitions would have important 
implications. In particular, one may combine this with information about the high tempera- 
ture phase transition for vanishing baryon density: One expects [15, 12, 14] an endpoint of the 
first order critical line between quark and nuclear matter if the zero density, high temperature 
transition is a crossover (as for two flavor QCD with non-vanishing quark masses). Such an 
endpoint corresponds to a second order transition where a large correlation length may lead 

1 The results of [12] indicate that condensates of quark Cooper pairs do not influence the behavior of the 
chiral condensate to a good approximation. The phenomenon of high density color superconductivity has only 
minor influence on the equation of state for quark matter [9] . On the other hand, for a substantial color octet 
condensate in the vacuum [10, 11] our assumption would not hold. 

2 Isospin violation and electromagnetism are important for nuclear matter in neutron stars. Our result for the 
equation of state is therefore not quantitatively realistic in all respects. Isospin violation can be incorporated in 
our formalism without conceptual difficulties. We have already included electromagnetism phenomcnologically 
for the quantitative description of nuclei. 
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to distinctive signatures in relativistic heavy ion collisions [16]. If the zero density transition 
turns out to be of first order in three flavor QCD, such an endpoint does not necessarily occur. 
(Endpoints are not excluded in this case, however, since the two first order regions could be 
disconnected.) The first order line for the gas-liquid nuclear transition exhibits a critical end- 
point for a temperature of about 10 MeV. Signatures and critical properties of this point have 
been studied through measurements of the yields of nuclear fragments in low energy heavy ion 
collisions [17, 18]. 

We emphasize that our treatment of quarks and baryons is still very crude and a different 
picture of the high density transition, e.g. due to the inclusion of other condensates, may well 
appear for two flavor QCD. A general outcome of our analysis concerns the crucial importance 
of confinement for any understanding of the phase transitions at high density and low temper- 
ature. In fact, the contribution of a free gas of baryons to the dependence of the free energy on 
the chemical potential is enhanced by a factor 27 (!) as compared to the contribution from a 
free gas of quarks. Therefore, the binding of quarks into baryons at low density and tempera- 
ture plays a crucial quantitative role which cannot be neglected by any satisfactory treatment 
of the high density transition. In fact, a first order transition between nuclear matter and quark 
matter would presumably connect an approximately free nucleon gas at low density to a quark 
gas at high density. Such a transition involves then a description in terms of baryons in the 
low density phase and cannot be understood within quark descriptions, which do not reflect 
the large baryonic enhancement factor. This may explain the phenomenologically unacceptable 
low critical densities often found in such quark descriptions. 

In quantum field theory the effects of a non-vanishing baryon density in thermal equilibrium 
or the vacuum are described by adding to the classical action a term proportional to the 
chemical potential /i, 

A^S = 3tfi £ b 3 j l J T dx J d 3 x ^7°^- = -3^5 . (1) 

The index j labels all fermionic degrees of freedom which carry a non-vanishing baryon number 
bj and a summation over spinor indices is assumed implicitly. For a description of the fermionic 
degrees of freedom in terms of quarks the sum is over N c colors and Np flavors, with bj = 1/3. 
We neglect the heavy quarks and concentrate on a two-flavor approximation where also the 
strange quark is omitted. For the nucleon degrees of freedom we include protons and neutrons 
with bj = 1. For our conventions, fi corresponds to the chemical potential of quark number 
density. The baryon number density n can be obtained from the //-dependence of the Euclidean 
effective action T, evaluated at its minimum for fixed temperature T and volume 3 V 



(B) i d r min T 

n = 



V 3dfi V 



(2) 

T,V 



We note that the Helmholtz free energy is F = r min T + 3fxnV. Our aim is a computation of 
the difference of r min between non-vanishing and vanishing /i. For T = this is dominated by 
fermionic fluctuations with (spatial) momenta q 2 < fj 2 . For not too large n (say /i ^ 600 MeV) 
we can therefore work with an effective model for the low momentum degrees of freedom of 

3 More precisely, B counts the number of baryons minus antibaryons. For T —* the factor T/V is simply 
the inverse volume of four-dimensional Euclidean space. 
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QCD. This argument generalizes to moderate temperatures, say T ^ 200 MeV. In the bosonic 
sector we will work with a linear meson model whereas for the fermions we keep the multiplet 
with lowest mass as discussed above. Our description takes into account the lightest scalar and 
pseudoscalar mesons as well as the lowest multiplet of vector mesons. 

The minimum of the effective action corresponds to the minimum of the effective meson 
potential U = TT/V for constant scalar meson fields. In consequence, U is a function of a 
complex Np x N-p scalar field matrix $, which describes the nonets of scalar and pseudoscalar 
mesons, and a similar matrix for the vector mesons. For a discussion of the chiral phase 
transition it will be sufficient to know the dependence of U on space and time independent 
fields which can acquire a vacuum expectation value consistent with SU(N F ) symmetry. These 
are the real diagonal elements of $ which we denote by a, and similar diagonal elements uj for 
the zero component of the vector mesons 4 . In the limit of vanishing current quark masses the 
minimum of U at sufficiently high temperature or high density should occur at a = in this 
model. For low T and \i spontaneous chiral symmetry breaking is triggered by a non-vanishing 
expectation value T), corresponding to the location of the minimum of U(a,u; //, T). (We 
adopt the convention through this work that bars indicate locations of potential minima.) The 
explicit breaking of chiral symmetry through non-vanishing current quark masses is described 
by a linear source term contained in U which induces nonzero a even in the phase without 
spontaneous symmetry breaking 5 [19, 21, 22]. The baryon density n, energy density e and 
pressure p follow from U (//, T) = U (<x(/i, T) , uJ(fi, T) ; /i, T) = e — Ts — 3 fin as 

n = -^U(fi,T); p — —U(fx, T) 

e = - = U(KT) + Zvn-T—{p,T). (3) 

Here we have normalized U (0, 0) = corresponding to vanishing pressure in the vacuum. 

For fluctuations in the momentum range < q 2 < (600 MeV) 2 we work within the linear 
quark meson model [21, 22, 23, 24, 25, 26] in an approximation which does not describe the 
effects of confinement. For low momenta, i.e. q 2 < q^, this description therefore becomes 
inappropriate. Three quarks are bound into color singlet nucleons. In this momentum range 
we describe the fermionic degrees of freedom by baryons, while keeping the description of the 
bosons in terms of the scalar field $ and corresponding vector meson fields. The use of the 
same bosonic fields for the whole momentum range will turn out to be an important advantage 
since it facilitates the computation of the free energy in different ranges of fi, corresponding in 
turn to different baryon densities and a different picture for the relevant fermionic degrees of 
freedom. For nuclear matter a typical value of the "transition momentum" is gn <^ 260 MeV. 
We find that the quark-hadron phase transition is substantially influenced by the change from 
quark to baryon fields at gn- This implies that a reliable quantitative understanding of this 
transition requires also a quantitative treatment of the change of effective fermionic degrees of 
freedom. 

4 The zero component of the vector fields lo^ can aquire a nonvanishing expectation value since at nonzero 
chemical potential Lorentz invariance is broken. 

5 For general quark or nucleon masses the diagonal elements a of the scalar meson matrix, and similarly 
for u>, can differ from each other. We suppress this dependence in the notation since we will only perform 
calculations for the case where all diagonal elements are the same. 
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Furthermore, our computation reveals that the transition from a nucleon gas to nuclear 
matter can be described realistically in terms of nucleon and meson degrees of freedom only if 
one accepts a relatively complicated form of the vacuum effective potential for the color-singlet 
chiral order parameter a. In particular, one needs large higher order couplings which do not 
seem very natural. It is conceivable that this situation changes for a more complex vacuum 
with additional condensates. 



2 Chemical potential flow equation 

We employ a new method for the computation of the /i-dependent part of the effective action 
that relies on an exact functional differential equation for V. This equation expresses the \i- 
derivative of T in terms of the exact field dependent fermion propagator. We start from the 
generating functional of the connected Green functions 

W[j] = In/ D X ex P j-S[ X ] - A,S[ X ] + / jx) (4) 

where x stands collectively for bosonic and fermionic fields with associated sources j and S 
is the action for /i = 0. For our purpose it is convenient to subtract from the effective action 
(defined by a Legendre transform) the /z-dependent fermion bilinear (1): 

r 5W 
r[p] = -Wkj\+J w-ApSfr], v> = -^-- (5) 

The /i-dependence of V arises only through A^S and can be expressed by a trace over the con- 
nected two-point function. Using the manipulations of generating functions outlined in [27] in a 
context with fermions [28, 29, 30] one obtains the exact nonperturbative functional differential 



equation 6 




dR * + rJ- 1 ) (6) 



d/j, 
where 

R»,jf(<l, q') = 3i/i6 i 7°(2vr) 4 5(g - q')5 jjf . (7) 

We remind that T is a functional of the meson and fermion fields, and the /z-derivative on the 
left hand side of (6) is taken for fixed fields. The exact inverse propagator is the second 
functional derivative with respect to the fields. It is a matrix in the space of internal indices 
and momenta and involves fermions and bosons. Since A^S* only affects fermions, the trace 
is over fermionic indices only and contains a momentum integration. For a configuration with 
constant bosonic fields and vanishing fermion fields does not mix bosons and fermions and 
is diagonal in momentum space. We therefore only need the inverse fermion propagator 

T^q') = H 3r m^)%q-q) (8) 

in order to obtain an exact equation for the /^-dependence of the effective potential 

9 Z = -£3&; / ^trry" \H( q ) + 3^/1 ~' . (9) 



dpi ^ J J (2 



in 



:i j 



3 We mention that in the presence of a local gauge symmetry this equation is manifestly gauge invariant. 



4 



Here tr denotes a Lorentz trace and Hjji, bjji = bjSjji are matrices in the space of fermion 
species. This exact relation expresses the baryon density for arbitrary quark mass term 7 (cf. 
eq. (30)) in terms of the exact fermionic propagators in presence of nonvanishing meson fields. 
We will see that the momentum integral is both ultraviolet and infrared finite such that eq. 
(9) is well defined. 

The exact fermion propagators are not known, and we have to proceed to approximations. 
The advantage of our underlying exact expression remains, however, that is easy to study 
which are the effects of qualitative and quantitative changes in the approximations for the 
fermionic propagators. In particular, we will learn how the transition from quark to nucleon 
degrees of freedom strongly affects the form of the effective potential - a discussion that would 
not be possible within a mean field approximation for a given effective model either of quarks 
or of baryons alone. In the present work we will use a rather simple approximation both for 
the quark and baryon propagators. For arbitrary a and uj we approximate 



We use our ansatz for the fermion propagator only to compute the /x-dependent contributions 
to the effective potential, i.e. we consider here the difference U(a, uj; /i, T) — U{a, uj; 0, T). In 
fact, the computation of the contributions due to a non-vanishing chemical potential allows 
one to use quite crude approximations in many situations. This is based on the observation 
that strongly interacting fermions are often successfully described as freely propagating quasi- 
particles. In our case they acquire an effective "constituent" mass m7 5 through a strong Yukawa 
coupling to mesonic vacuum expectation values. (The matrix 75 appears in the mass term 
as a consequence of our Euclidean conventions [28].) Similarly the constant field uj denotes 
the analytic continuation of the zero component of the Euclidean u>-vector-meson field with 
coupling to the fermions. The piece ~ Q is the remnant of the vector coupling ~ uj^ in a 
situation where the zero component of uj^ can acquire an expectation value due to the breaking 
of Lorentz invariance 8 by the nonvanishing chemical potential (1). We will later determine the 
values of uj self-consistently. An important simplification in our ansatz (10) is the neglection 
of a possible wave function renormalization Z^(q,a,uj) which could multiply the kinetic term 
q u Y '. In a more realistic setting this will certainly play a role near the transition between 
quarks and baryons. 

We emphasize that the computation of dU/d/i according to eq. (9) does not need any 
information about the masses and effective self-interactions of the mesons. They determine, 
however, the effective potential at zero baryon density U(a, u;0,T) and therefore influence 
the possible phase transitions. In fact, the meson self-interactions may turn out to be quite 
complicated. We do not attempt here to compute the meson masses and self-interactions by a 
mean field approximation, since earlier renormalization group investigations have shown [22] 
that this is probably much too crude. The advantage of our method is that the lack of knowledge 
about the meson interactions can be separated from computation of the /i-dependence of the 

7 The relation (9) holds for arbitrary a which corresponds to arbitrary quark mass m q through dU/da ~ m q . 
8 For fi — one therefore has ui = 0. 



Hjj'(q) = q u "f u + mj(a; fi, T)>y 5 + ibjVLrf 5 jjt 



(10) 



with 




(11) 
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potential. We will simply parametrize U(a, u; 0, T) in accordance with the symmetries and 
(indirect) observational knowledge. 

At this point it may be useful to summarize the approximations that affect the computation 
of the /^-dependence of the effective potential if we use the ansatz (10) in the exact relation 
(9). Perhaps most importantly we omit the dependence of the fermion wave function renormal- 
ization on momentum, a, u, fi and T. We also neglect a possible difference in normalization 
of the quark kinetic term and the baryon number current. Similarly, we have not considered 
a possible momentum dependence of the mass term as well as the momentum dependence of 
the contribution ~ 7 . Finally, we assume that m,j can be taken as independent of u, and g u 
as not dependent on uj or ji. In this approximation the term ~ Q can be combined with 
such that ji is replaced in the propagator (r( 2 ) + R^ 1 by an effective chemical potential 



yUeflf = A*+ -fi(cj,cr;T) 



(12) 



With the approximation (10) the evolution equation for the /x-dependence of the effective 
meson potential takes a very simple form 



dU 
d/j, 



(2k) 



The remaining trace over spinor indices is easily performed 

dU 
d/j, 



— — = — 2 



K 3 = 



E / d3 ' 

j 

3 J 2tt 



K 

(2tt)3 3 ' 



(q° + Sibjfies) 2 + q 2 + m 2 j 



(13) 



(14) 



(15) 



For non-vanishing temperature the g°-integration is replaced by a sum over Matsubara fre- 
quencies 

dq° 



^ >TY 

^ ne2Z 



(16) 



with q° = 2n(n + 1/2)T and, correspondingly, 5(q — q') — > 5(q — q')5 nn >/(2nT). This results in 



+ 1 



e T + l 



(17) 



where the two terms are proportional to the fermion and anti-fermion contributions to n. We 
see explicitly that the momentum integration is finite due to the exponential suppression for 
large if. 

We will concentrate here mainly on T = where the g°-integration yields a step function: 



Kj = 3&,e(9& 2 (M 2 - if + m))) . 



(18) 



The remaining g-integration is therefore cut off in the ultraviolet, q 2 < 9b 2 (fi c g) 2 — m 2 , and 
only involves momenta smaller than the Fermi energy 36^^. As it should be, it is dominated 
by modes with energy Jq 2 + m 2 near the Fermi surface. 
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3 Quark and nucleon degrees of freedom 



We will assume that dU/d/j, can be expressed as a simple sum of the contribution from quarks 
with momenta q 2 > <2h and that of baryons with momenta q 2 < q^. This is the simplest 
approximation which catches the effective transition from quarks to baryons as effective degrees 
of freedom in the relevant momentum range. It will be sufficient to demonstrate the most 
important effects of confinement on the /i-dependence of the effective potential, namely that 
the contribution of a gas of nucleons is greatly enhanced as compared to a corresponding 
contribution of quarks. In a more realistic scenario the transition between quark and nucleon 
degrees of freedom will be less abrupt. Within our approximations part of the uncertainty 
related to this effective transition can be studied by allowing that qn depends on a since 
typically the relevant values for a depend on the baryon density, being higher for a nucleon 
gas than for a quark gas. 

We first consider the range of momenta with q 2 > q^ for which we use an effective linear 
quark meson model [21, 22, 23, 24, 25, 26]. Here the quark mass term ~ 7 5 arises through 
a Yukawa coupling h to the expectation value of the cr-field, m q = h q (a; fi, T)a. Since the 
quark description breaks down for small momenta, we restrict the integration over q 2 in (14) 
to the range q 2 > q^. We therefore infer for the quark contribution to the /z-dependence of 
the effective potential (for /i > and T = 0) 



9U (Q) N C N F I"/ , , 2 2\3/2 -2 2 „2 \ (, Q S 



For the low momentum range q 2 < q^ where the fermionic degrees of freedom are the 
lightest baryons rather than quarks we repeat the steps leading from (6) to (19). The trace 
now involves a sum over proton and neutron but no color factor. This yields a contribution 
(again for T = 0) 



0l/(B) 



(9/4 - m^) 3 / 2 0(9^ ff _ m 2 N )Q{m 2 N + q 2 H - 9/i e 2 ff ) 



djJL TV 2 ^ 

+ q 3 H Q(9fi 2 cfi -m 2 N -q 2 H )y (20) 

We parametrize the nucleon mass as m^cr) = 3/ijv(cr)<7 and note that for ~ h q one has 
mjy — 3m g , as appropriate for nucleons described as composites of three constituent quarks. 
The baryon density can be directly inferred from eqs. (19), (20) as 



_i f duw du^ \ 

3 \ dfj, d/j, J 



(21) 



cr=cr, W=LU 



since the partial derivatives of the effective potential with respect to a and uj vanish at the 
yU-dependent potential minimum (a, a;). 

We repeat that this picture is only a crude approximation to the binding of quarks into 
nucleons. A nucleon description should work well for h^a 2 near nl s , since only low momentum 
degrees of freedom contribute in this range. On the other hand, the quark description becomes 
important for h 2 a 2 n 2 s — q^. In a more realistic description the 6-functions in (19), (20) 
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would become smooth. The characteristic quark-baryon transition momentum gn typically 
depends on a. Indeed, a baryon description for the low momentum degrees of freedom is 
necessary for a not too far from its vacuum expectation value a . We will see that in this 
range of a the transition momentum q H is around 260 MeV or higher. On the other hand, 
baryons do not seem to be meaningful degrees of freedom in a situation of chiral symmetry 
restoration at a = such that qn may vanish for a = 0. 

For hj<s of the same order as h, the nucleon mass is about three times the quark mass for 
a given value of a. Therefore eq. (20) results in an important enhancement of dU/dfi in the 
range ^ fi 2 s — < fo N \ a\ < fi c s as compared to the contribution from the quarks. This 
is mainly due to the fact that more energy levels fall below the Fermi energy 3/z e fj for the 
baryons. More precisely, a factor 3 3 = 27 arises from the ratio [(V - m 2 N )/(^ 2 - h 2 q a 2 )f/ 2 
if the subtraction of q z H in eq. (19) can be neglected. An additional suppression of the quark 
contribution by a factor b q = 1/3 in the coupling of the chemical potential is canceled by the 
color factor N c = 3. (Neglecting strangeness, the two quark flavors are matched by the two 
species of nucleons. For a light strange quark one would observe an additional enhancement 
for the baryon contribution due to the larger number of baryons in an octet as compared to 
the three flavors.) This "nucleon enhancement" is one of the most important observations of 
the present paper. We believe that this effect is quite robust in view of a possible more precise 
modeling since only very simple properties of the fermion degrees of freedom play a role for our 
argument. In our description the large "nucleon enhancement" by a factor of about 27 is the 
basic mechanism which may lead to separate gas-liquid and hadron-quark phase transitions. 
Despite this enhancement one observes that dU /dfi is continuous in a and fi. Furthermore, for 
^at(c) = h q (a) the simultaneous jump of the renormalized fermion mass by a factor of three, 
together with a similar jump of the renormalized coupling to vector mesons (due to the factor 
bj in eq. (10)) could also be accounted for by a sudden drop of the fermion wave function 
renormalization 9 form = 1 for q 2 > q 2 H to Z^ = 1/3 for q 2 < q 2 H . This corresponds to the 
continuity in /x eff which does not depend on the wave function renormalization multiplying the 
fermion kinetic term. 

With qj^- = JL we can easily rewrite eqs. (19), (20) as flow equations for /i c g. In the 
approximation of /i-independent Yukawa couplings h = h(a), = h^(cr) and gn = <?h(c) 
these differential equations can be integrated analytically. We define 

U{a,w;n,T) = U (a;T) + U u (a,u;T) + 2U li (a,u; f jL,T) (22) 

where 2U^(a, uj; fx, T) entails the /z-dependent contribution from the two lightest quarks (2U^) 
as well as proton and neutron (2U^). The /^-independent part of the potential is thus given 
by Uq + U u with U u the u;-dependent contribution. For T = we obtain 

U, = UF + Up), (23) 
/^cflf (rfe - \ h2 q a2 ) vVcfr - h 2 q ° 2 

9 Such a drop of the effective wave function renormalization would be required for a Higgs picture of the 
QCD vacuum where quarks and baryons are described by the same field [10]. 



E^(<7,c^,0) = 
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+ 



3,4 4 , ^ff+y^ff-^ 2 




(24) 



g H (4g^ eff - (3g 2 + lhy)Jh*o* + &) 1 6(/i c 2 ff - ^a 2 - g 2 ) , 



UV>(a,u>;iJL,0) 




Vesiffis ~ -/iN^ 2 )\/^cfr - h l° 2 



+ 




+ 



+ 




(25) 



In this expression the dependence on /x and cj appears only implicitly through /i cff . For given 
h q and /i^r the a-dependence of is uniquely determined once the dependence of qn on a is 
fixed. 

The qualitative dependence of qu(cr) on a can be inferred from the following argument: A 
crucial ingredient for the confinement of quarks in hadrons is the formation of QCD strings. 
Strings break because of pair production of mesons if typical quark kinetic energies become too 
large. Therefore baryons can only exist for sufficiently small average quark kinetic energies or 
momenta. Very roughly, the relevant critical kinetic energy is expected to be proportional to 



the pion mass y <?h(<7) + h 2 a 2 ~ 2m 7T (a). The a-dependence of the pion mass can be inferred 
from the effective potential as m 2 (a") = (dU/da + 2m 2 / 7r )/(4cr), with m n the pion mass in the 
vacuum and f n the pion decay constant. Since m 2 (cx) always tends to zero for small enough a 
there should be a critical value a c for which gn( cr c) = 0. We will not use baryons for a < a c 
and take qn(a < a c ) = 0. For our purpose we will be satisfied with a crude approximation 10 
where we neglect the a-dependence of gn m the range of a relevant for nuclear physics, a > o" H 



quia) = -^^yV 2 - °l) (2^1 - ol - a 2 )Q(a - a c )0(a H - a) + q u Q(a - a u ) . (26) 



For a wide range of a c and <th our results for the nuclear gas-liquid transition will turn out to be 
independent of the precise values of these two quantities. For definiteness we take a c — 15 MeV, 
o"h = 25 MeV. We expect that the constant qn should have the size of a typical QCD scale, i.e., 
around 200 MeV. On the other hand, we will find that the quantitative aspects of the quark- 
hadron transition depend on a c , o"h and gn which parameterize in our crude approximation the 
effects of confinement. This underlines that a more quantitative understanding of the effective 
transition from quarks to nucleons is needed before reliable statements about the hadron-quark 
phase transition at low temperature can be made. 

It is interesting to note that for two light quark flavors (N c Np = 6) the /i-dependent 
contribution to the potential at the origin and therefore to the energy density reads for arbitrary 

10 For our choice dqn(a)/da is continuous at a = <7h- 
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finite h(a, //) 

ef = ~QU,(a = 0,^,0) = ^a4 = > (n^f* . (27) 

This has the simple interpretation of the total energy of six massless quarks with all energy 
levels filled up to the Fermi energy /x e g. Furthermore, for a and /i in the range relevant for 
nuclear physics and for sufficiently large g H , i-e. > 9(nls ~ ^<J 2 ), the contribution Uj^ is 
simply the mean field result for a nucleon meson model, whereas Uffi vanishes. Our approach 
gives a new motivation for the approximate validity of mean field theory from the truncation 
of an exact flow equation. Furthermore, it offers the possibility of a systematic improvement, 
e.g., by taking the /i-dependence of into account. Despite this similarity, our method goes 
beyond mean field theory in an important aspect: For the free energy only the difference 
between vanishing and non-vanishing chemical potential is described by mean field theory, 
whereas we do not rely on mean field results for the effective action r at \x = 0. Since T(/i = 0) 
is the generating functional for the propagators and vertices in vacuum it can, in principle, be 
directly related to measured properties like meson masses and decays. This is very important 
in practice, since mean field theory does not give a very reliable description of the vacuum 
properties. 

4 Meson interactions 

In order to discuss possible phase transitions as fi is increased beyond a critical value we 
need information about the effective potential for /i = 0. For a vacuum without spontaneous 
symmetry breaking relatively accurate information about Uo(a; T) — Uo(a, uj — 0; /i — 0, T) for 
all relevant a could be extracted from the knowledge of meson masses and interactions. Also the 
approximation (10) for the fermionic propagator would presumably be reasonable for arbitrary 
a. In case of spontaneous chiral symmetry breaking the situation is more complex: The true 
effective potential U becomes convex because of fluctuations which interpolate between the 
minima of the "perturbative" or "coarse grained" potential [31, 32]. Masses and interactions 
give only information about the "outer region" of the potential which is not affected by this 
type of fluctuations. In parallel, the simple form of the fermionic propagator (10) becomes 
invalid in the "inner region" for small a because of a complex momentum dependence [31, 32] 
and the breakdown of the approximation of a constant Yukawa coupling. In order to cope with 
these difficulties, U (a; T) should rather be associated with a coarse grained effective potential. 
For a suitable coarse graining scale 11 k the effect of the omitted fluctuations with momenta 
smaller than k is expected to be small near the /i-dependent minimum of U. Around the 
minimum at a we can therefore continue to associate Uo(a;T) with the effective potential 
and relate its properties to the measured masses and decay constants. On the other hand, we 
do not have much information about the shape of U (cr; T) for a ~ 0. This uncertainty in the 
appropriate choice of Uo(a; T) is one of the main shortcomings of our method. In practice, we 
interpolate the partly known polynomial form of Uo(a; T) form the outer region (which includes 
the minimum characterizing the vacuum) to the inner region for small a. By continuity, this 

n The coarse graining scale k is chosen such that is approximately fc-independent for \<j\ around \a\ or 
larger, whereas the approach to convexity for |<r| < \a\ and k — > has not yet set in. 
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should be quite reasonable for nuclear matter since the relevant values of a are not much 
smaller than the vacuum expectation value a . For quark matter, the uncertainties are more 
important. 

We investigate here the two flavor case with a potential of the form 



U (a;T) 



+ - 



4 7 3 (T) 
3<7 2 (T) 
4 7s(T) 
5a§(T) 



a 2 — (Jq(T) +2A(T) 



a 



74(T) 



"'SCO 



+ 

4 



n5 



^-^(T) -2 J( t + c(j,T) 



where 



j = 2ml(0)a (0), c(j, 0) = 2ja (0) 



(28) 



(29) 



In the remainder of this work we mainly consider T = and use A = A(0), U (er; h) = U (a; fi, 0) 
etc. The meson field is normalized such that a = o"o(0) is related to the pion decay constant 
by (T = /tt/2 = 46.5 MeV. This means that the pions have a standard kinetic term (as derived 
from Cfo n = Tr d^d^Q). Because of higher order kinetic invariants [24] the kinetic term 
for the sigma meson, £kin,o- = 2Z a d^crd^a can involve a wave function renormalization Z a 
different from the one for the pions. The potential (28) arises from a fifth order polynomial 
in the invariant p = Tr$^$ = 2a 2 with an additional source term —^Trj(<& + where j is 
proportional to the renormalized current quark mass (say at 1 GeV) . The only violation of the 
chiral SUl{2) x SUr(2) symmetry arises from this source and in the chiral limit of vanishing 
current quark masses the last two terms in eq. (28) should be dropped. The coupling A is 
related to the a-mass m a by fh 2 a = Z a m 2 = m 2 + 4A<Tg. We will use here m a = 510 MeV, 
A = 28. It is actually m CT rather than the physical mass m a which is relevant for the properties 
of nuclear matter. One of the parameters 73, 74 or 75 can be eliminated in favor of the scale 
/Iq which characterizes the hight of U at the origin 



t4 - 



2tt 2 



= 17(0; 0) = 4 (2X - -73 + 74 - - 75 ) + 



4^0 • 



(30) 



Without the complications of confinement (i.e., for gn = 0) the quark-hadron phase transition 
in the chiral limit (j = 0) would occur for /i e s = /x . 

Finally, we determine the expectation value of ui by observing the identity 



du 



2 dU„ 



dfx 



(31) 



It follows from first differentiating eqs. (19), (20) with respect to uo and then performing the 
/^-integration. For the /i — contribution we only take into account a a- and T-dependent 
mass term 

(32) 



U^ = --M 2 (a,T)iu 2 . 



The solution of the cj-field equations for arbitrary a, T, /i obeys 

u(a, fi, T) 



l9w 9Uf, 



(33) 
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We note that at the potential minimum oJ is proportional to the baryon density with a negative 
coefficient. This implies that the coupling to u reduces the effective chemical potential [13]. 
In the following we will always assume that oJ(cr, fi, T) is inserted such that fi e s becomes a 
function of a, /i and T. The field equation which determines the location a of the minimum of 
the effective potential can be expressed in terms of partial cr-derivatives at fixed /z e g- 

+ W=0. (34, 

Below we will also neglect a possible cr-dependence of M w . The location of the minimum a 
becomes then independent of the value of uJ and only depends on /x e g. For this setting the 
coupling to vector mesons is relevant only for the relation between /x eff and /i. 



5 Meson— baryon interactions 

A crucial ingredient for any quantitative analysis is the sigma-nucleon coupling h-^- We first 
investigate if chiral symmetry and the observed value of the pion nucleon coupling place any 
restrictions on this coupling. For this purpose we employ a derivative expansion of the most 
general effective Lagrangian which is bilinear in the nucleon doublet field \1/n and involves 
scalar and pseudoscalar fields contained in the 2x2 matrix $ 

L = ^|^NRi 71 ($$ t ,p)$^NL -"^N L i 71 ($ t $,p)$ t ^NR 

+ * NL G'i($$ t ,p)i7^ /1 *NL + *NRG'i($ t $,p)i7^NR (35) 

+ ^NL^G^^p)^^ • 

Here we have imposed V and C symmetry and used \I/ NL = (1 + 7 5 )^n/2. With the standard 
decomposition 

$ = a!; 2 = all , £ = exp ( -7-7*7? 



one finds 



Aa J (36) 



C = 3h N (a)aNi 5 N + Z N (a)N(i^dr-^Vp + G A ((TW'fa li )N 
.Z N (a) 



2a 2 



[G A (a) - l]a(d lx a)N^N (37) 



where 



h N (a) = F(a 2 ,2a 2 )/3 , Z N (a) = G x {a 2 ,2a 2 ) 

2G 2 (a 2 ,2a 2 )a 2 (38) 



G A (a) = 1 



Gi(a 2 ,2a 2 ) 
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and 

= ~{^ + ^) 
a, = - l -(^d^-id^) = ±-Td^+.... (39) 

Normalization of the baryon number current requires Zn(o"o) = 1 and we neglect the a- 
dependence of in the following. The strength of the linear pion-nucleon coupling is fixed 
by #a = Ca(o"o) an d bares no relation to the function h^(a). We may expand h^(a) around 

h N (a) = h N (a ) + ^ (a 2 - a 2 ) + . . . . (40) 
°o 

With t«n = 3/iN(c"o)o"o = 939 MeV we find h^(cr ) = 6.73. Linearizing m^(cr) = 3/in(o"V 
around <7 then yields m N (<r) = 3hcr + e G with h = h^(a ) + 2g N , e& = —6g^cr . The linear 
sigma-nucleon coupling h is a free parameter which is expected to be in the vicinity of h N (a ). 
We will determine it below from the properties of nuclear matter. Since h also appears in the 
scattering of nucleons a comparison with experiment may serve as a test for our model. 

6 Nuclear matter and the nuclear phase transition 

Let us turn to the zero temperature properties of nuclear matter in our picture. For ji = the 
effective potential or free energy U has its minimum at a = f w /2 = 46.5 MeV. The potential 
in the region near o"o is not altered as long as // e ff remains small enough (cf. eq. (20)). This 
changes as fi c s is increased beyond a critical threshold. For suitable parameters in Uq (eq. 
(28)) we observe that for 3/i somewhat below the nucleon mass a new minimum of U occurs 
at a( nnc \fi) < a , with a potential barrier between both minima. For a certain range of /x the 
local minimum at a^ nnc \fi) and the global minimum at a coexist. As /x increases, the value 
of {7((f( nuc )(/x)) is lowered whereas U(a ) = remains fixed as long as the effective chemical 
potential is smaller than a third of the nucleon mass, /x e g < /in(c"o)°"o- There is a critical value 
/x nuc for which the two minima at 

a nuc = a^ c \fi nuc ) (41) 

and ctq are degenerate, U(a nuc , /x nU c) = ^(co,A*nuc) = 0. The corresponding critical potential is 
plotted in figure 1. Both phases have equal, vanishing pressure p = —U and can coexist. We 
observe that the phase transition between the vacuum (a = <t ) and nuclear matter (a = (x n uc) 
is clearly of first order. For small temperature this corresponds to the transition between a gas 
of nucleons and nuclear matter which may be associated with a nuclear liquid. 



For a quantitative description we concentrate mainly on large values of qu where this 
transition happens in a region with /ig ff < ^N^nuc^nuc + <?h/9- In this case the transition 
from nucleon to quark degrees of freedom does not affect nuclear properties and the phase 
transition from a nucleon gas to nuclear matter. We partly recover the cr-u;-model of nuclear 
physics [33, 34, 20], in a context where chiral symmetry breaking and constraints from meson 
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35 40 

a I MeV 

Figure 1: The critical effective potential for the nuclear liquid-gas transition corresponding to the parameter 
set A of table 1. 



masses and decays are properly incorporated. For large enough qu the critical baryon density 
of the nuclear liquid is given by eq. (20) 



7l I: 
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71- 



/^cff^nuc) °"nuc) ~~ h N ((7 DUC )(J 



N T 2 

nuc 



3/2 



(42) 



We will see below that one can identify n nuc with the baryon density in nuclei = 
1.175 x 10 6 ■ MeV 3 up to small corrections. Furthermore, the baryon number indepen- 
dent contribution to the binding energy per nucleon in a large sample of nuclear matter 
is known from the mass formula for nuclei: (3^ = —16.3 MeV. In our context one finds 
(3^ = 3/i nuc — ttin and for realistic models the gas-liquid transition should therefore occur 
for yUnuc = 307.57 MeV. Eq. (42) then yields a quantitative relation between the effective 
chemical potential in nuclear matter // e fr,nuc = A*eff(A f nuc, <?nuc) an d the effective nucleon mass 



3/i N (cT r 



<Jr 



€q + 3ha nuc . For m N (a nuc )/m^ = (0.6,0.65,0.7,0.75,0.8) one 



m N (a nuc ) 

finds /icfr,nuc = (206.6, 220.9, 235.4, 250.1, 264.8) MeV. Equivalently, this can be seen as a rela- 
tion between h^(a nuc )a nuc and the coupling if we use with (34) 



3M2' 



(43) 



For a (T-independent u;-mass M u = 783 MeV typical values for the above ratios for 
m^(a nuc )/m^ are = (12.61,11.68,10.66,9.52,8.21). From the value of nuclear density 
we can compute the Fermi momentum g nuc = 259 MeV. This yields for this scenario a lower 
bound gn > <?nuc = 259 MeV. (For quantitative computations we take gn = l-2g n uc-) 
An important quantity for the equation of state is the compression modulus 



K 



9n z 



dri 2 



n 



9(^-2?- 
\ an 



n , 



(44) 
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h 




73 


74 


[5 


Mo 

McV 


A 


5.4 


9.02 


-30 


47 


-60 


372 


B 


5.0 


9.52 


19 


112 





348 


C 


4.6 


8.74 





55 





330 



Table 1: Coupling constants for three different parameter sets. The linear sigma-nucleon coupling h, the 
coupling of the w-meson to the it, d-quarks and the nuclcons, the meson self-interactions 73, 74, 75 and 
the scale /zq are defined in sections 4 and 5. 



The value at the phase transition 
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9^ 
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+ 9n n 



^nuc 

(w)2 



97 ^ 
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M 2 

0J 



+ 



^nuc 



2/3 



(inn 



hm^(a nuc ) 

/^eff,nuc 



da 



(45) 



has been inferred from experiment as K = (210 — 220) MeV [35, 20]. This can be used to obtain 
additional information about the //-independent part Uq of the effective potential. Neglecting 
the a-dependence of M w one finds by differentiating eq. (34) 



da 
dn 



hm N (a) 



~da^^ 



+ 2 

da 2 



6 ~ 2 2 V^eff-^N^) 



(a) H -/i m N (a 



^eff 



7T Z 



/^eflf 



(46) 



in 



Combining eqs. (46) and (45) the compression modulus yields information about ^g^-(a n 
addition to Uo(a nuc ) and ^-(a nuc ) which are determined (for given mN^nuc) and h) by the 
condition U(a) = and the field equation (34). 

For any given value of the coupling 75 the system of equations provides a mapping between 
the parameters (h, 73, 74) and the quantities (n nuc , /3, K , m^(a nuc )). For a demonstration 
of the range of values for various quantities of interest we report our results for two parameter 
sets with different 75 (A and B) in tables 1-3. (For both sets (5 = —16.3 MeV and n nuc = n^ n \) 
Agreement with nuclear properties can indeed be achieved. It is not our aim here to make 
a precise determination of parameters and we only mention that somewhat smaller values 
of m N (a nuc ) or other (large) values of the compression modulus lead to qualitatively similar 
results. One finding remains common, however: for nuclear matter properties in a reasonable 
range we always need large values of some of the couplings | j , |7 4 | or |7 5 |. No viable solution 
was found for an approximately quartic meson potential with small 173,4,5). This may be a cause 
of worry for this class of models since earlier renormalization group studies of the effective 
meson potential have typically resulted in substantially smaller higher order couplings (73,4,5 1 
than the ones needed here [22]. 

For given parameters we can also compute the energy density and the pressure and relate 
it to the baryon density. This determination of the equation of state of dense nuclear matter is 
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rather insensitive to aspects of QCD not treated here like the role of gluons or vector mesons 
(beyond the effect of uJ 7^ 0). The reason is that these degrees of freedom do not contribute 
to the difference of the effective action between zero and nonzero chemical potential. The 
nuclear equation of state can therefore be considered as a prediction of the model (for fixed 
parameters). 





& nuc 

MeV 


m- N (o- nuc ) 
m N 


K 
MeV 


£ 


Mcff .nuc 

MeV 


A 


33.2 


0.77 


214 


1 


256 


B 


30.85 


0.75 


217 


1.3 


250 


C 


29.8 


0.755 


oo 


2 


259 



Table 2: Properties of nuclear matter at vanishing pressure for the parameters of table 1. The table shows 
values for the chiral order parameter cr nuc , the effective nucleon mass m^{a mc ), the compression modulus K n 
and the effective chemical potential /i ff.nuc for nuclear matter. The surface tension S for the droplet model of 
nuclei normalized to the value £( n ) extracted from the nuclear mass formula is discussed in section 7. 

In figure 2 we have plotted the binding energy per nucleon, f3 = e/n — m^, as a function 
of density corresponding to the parameter set A. For values of n larger than approximately 
1.7 the details of the transition from nuclear to quark degrees of freedom become important 
and we don't expect our results to remain quantitatively reliable. Similarly, our results for 
the baryon density as a function of pressure are displayed in figure 3. Figures 2 and 3 can be 
combined to yield the equation of state e{p) for n < 1.5ri^. 




1 1.1 1.2 1.3 1.4 1.5 

Figure 2: Binding energy per nucleon (3 as a function of density. Parameter values correspond to A in table 1. 
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Figure 3: Baryon density as a function of pressure in the vicinity of the nuclear gas-liquid transition at very 
low T. Parameters correspond to set A in table 1. 

Within the " a — uj model" (models A and B) the dominant repulsion between nucleons at 
short distance is ascribed to the exchange of a>mesons. It is not established if this repulsion is 
indeed sufficient. Another possible repulsion mechanism is the effective transition from nucleons 
to quarks at short distance. In view of the potential difficulties of the a — uj model we also 
explore this second alternative — our model C. For this purpose it is instructive to consider 
an extreme scenario where the characteristic quark-baryon transition momentum q H takes on 
its lower bound 



The results correspond to the set C of tables 1-3. Because of the G-function in eq. (20) the 
nucleon contribution to the density does not increase any more as /x e fr exceeds the critical 
value given by eq. (43). On the other hand, there is a range of /i e s for which the quark 
fluctuations (19) do not yet contribute to the baryon density. For this range the density n will 
not depend on any other parameter of the model and n nuc = is guaranteed by eq. (47). 
Details of the potential in the vicinity of a nuc are now affected by the transition from nucleon 
to quark degrees of freedom. For the simple choice (26), however, they do not depend on a c 
or <th provided both are smaller than <7 nuc . Because of the gap in fj, e g between the nucleon 
Fermi surface and the onset of quark fluctuations many properties become very simple. The 
minimum occurs within the range a q < a QUC < a n {. Here <r q corresponds to the onset of quark 
fluctuations 



whereas cr n f denotes the maximal value of a for which all nucleon levels with q 2 < are filled 




(47) 




(48) 




(49) 
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For values of a^ nuc \fx) between max(a q ,au) and <7 n f the baryon density is independent of fx 

n(fx) = ^qg = . (50) 

This would give a natural explanation for a large compression modulus K according to eq. 
(45). Also, for max(cr q , <r H ) < o < <7 n f one finds that dU^/dfx is independent of a and the 
constant shift (43) between fx e g(a, A 4 ) and /z holds for all a . The relation between n nuc and is 
such that up to a Fermi momentum (^(ow) = all levels are filled with nucleons (or bound 
quarks). In this crude picture the higher momentum levels (corresponding to a larger baryon 
number in a fixed volume) would have to be filled by free (constituent) quarks. This leads 
to a particularly simple explanation why nuclear density is almost independent of all other 
parameters characterizing the state of nuclear matter at T = 0, like pressure, baryon number 
or the Z/B ratio of a nucleus. Typical parameter values and corresponding characteristics of 
nuclear matter for this "saturation scenario" can be found as set C in tables 1 and 2. 

We next discuss the equation of state for the saturation scenario. For a given value of a 
one finds in this scenario a range fx n f < fx e s < fx q with constant dU^/dfx where 



M^) = y/h*o* + <&(*). (51) 
In this range uj is independent of fx and has the simple form 

U^a; fx) = U^a- ~p) ~2 (a* A*) <&{<?) • (52) 

Here fx is a fixed reference value within the interval [/x n f, /x q ], and we remind the reader that 
quic > ctr) = qu- We note that the location of the potential minimum at a^ nnc \fx) is indepen- 
dent of fx. Using ~jl = fx nuc , the pressure and energy density of nuclear matter are 

3 

p = 2%/i-/i nuc ) (53) 
2 

6 = -^Qnf 1 ~P = ( m N + P)n . (54) 
The (T = 0) equation of state for nuclear matter for this extreme saturation scenario 

dp dp ^ ' 

implies a diverging compression modulus Kq and is therefore not fully realistic. Nevertheless, it 
is well conceivable that the true behavior of nuclear matter is somewhere between the simple 
version of the a-cu model and the extreme saturation scenario. In the language of the a-uj 
model this would be expressed through the momentum-dependence of couplings and wave 
function renormalizations (form factors). 

Despite the substantial difference in the compression modulus the three scenarios (A)-(C) 
all show a similar value of <7 nuc ~ 30 MeV and therefore a nucleon mass in nuclear matter 
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around 700 MeV, as supported by some experimental evidence [36, 20]. Also the critical value 
/-teff.nuc — 250 MeV is very similar for these three models. All three scenarios support the 
existence of a first-order transition from the vacuum to nuclear matter. Extending this result 
to small T ^ this results in a first-order transition from a nucleon gas at low density to a 
nuclear liquid at high density, in complete analogy to the vapor-water transition. 



7 Droplet model for nuclei 

For a first-order gas-liquid transition we can describe sufficiently large nuclei as droplets of 
nuclear matter in a surrounding vacuum. For quantitative estimates of their properties we 
have to take into account that because of the surface tension the pressure inside the droplet 
is different from zero. The nucleus is at equilibrium if the pressure equals the derivative of the 
sum of surface and Coulomb energy 12 with respect to the volume 

dEx , dE c 



P dV + dV 

1 /3 

£ s = (SQ^/^V 2 / 3 , E c = ^-(^j kZV- 1 ' 3 . (56) 

The surface tension E can be expressed (thin wall approximation) in terms of the potential as 

Z \ 1/2 



<J nuc 



E = 2 -M / daJ2U(a-»). (57) 



We use the phenomenological relation Z/B ~ (2 + 0.01535 2 / 3 ) a — 1/137, and k should be 
very close to one. This leads to a volume- and therefore baryon number dependent pressure 

p= i(lf^B-*-±(lf aK <L*ir* (58, 

and a total binding energy per nucleon 

Eb n , x p Ex E c n £ s E c 9 p 2 

Neglecting the pressure term the mass formula for nuclei yields the "experimental" values 

E (n) = 4.22 -10 4 MeV 3 

n (n) = 0.96 . (60) 

For B = 208(12) one finds p = 0.9(5.9) • 10 6 MeV 4 . The pressure therefore contributes to E/B 
only very little, A(E/B) = — 0.012(— 0.53) MeV and can indeed be neglected for large B. For 
small B eqs. (58), (59) result in an interesting correction to the mass formula, which is usually 
not taken into account in the droplet model for nuclei. Comparison with fig. 3 shows that 
our model yields a baryon density which is indeed almost independent of B for large nuclei. 



12 



We neglect here the asymmetry effect from the proton-nucleon mass difference. 
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(Note that formally corresponds to B — > oo.) On the other hand, for small nuclei the 
baryon density is enhanced. The value of the surface tension as computed from eq. (57) for 
different parameter sets can be found in table 2. This result is quite reasonable in view of the 
uncertainties, first from the proper choice of a coarse grained potential in (57) (cf. [37, 38, 39]), 
and second from the choice of parameters in U . In fact, the successful explanation of the small 
ratio (S^) 1 / 3 / g nuc is encouraging. In summary, our simple approach gives a quite reasonable 
picture for nuclei. We emphasize that once Uo(a) is fixed, our approximations allow for a "first 
principle calculation" of properties of nuclei! 

8 Two flavor quark hadron phase transition 

At the critical chemical potential // nuc the free energy U (er; /i) shows two degenerate minima: 
one at <7o, with vanishing density and the other at <7 nuc , where nuclear matter density n nuc is 
reached. In this section we consider densities higher than n nuc . For sufficiently high density 
one may expect, and we observe, a further transition from nuclear matter to quark matter. A 
new first-order phase transition would be related to a third distinct minimum of the effective 
potential U. 

The results of a quantitative analysis for the polynomial potential (28) with the parameter 
sets A-C are reported in table 3. Using for the transition momentum the ansatz of eq. (26) 





A'qm 


/^eff ,qm 


Tlq m 


<Tq m 




y 


MeV 


MeV 


^nuc 


MeV 


MeV 


10 6 MeV 


A 


975.5 


530 


8.6 


2.2 


15.8 


5.9 


B 


859.7 


484 


6.5 


2.2 


15.4 


4.8 


C 


511 


370 


2.9 


5 


34 


1.7 



Table 3: Critical quantities for the quark-hadron phase transition. The values for the chemical potentials /^ qm 
and /^ e ff,qm, the baryon density in the quark matter phase n qm , the order parameter <7 qm , the effective quark 
mass M qiqm and the surface tension S qm for the sets A and B should be interpreted as an illustration of the 
uncertainties of a polynomial extrapolation of the potential Uq to the origin a = 0. 

with a c and er# in a reasonable range, we find a first order transition between nuclear and 
quark matter. The sets A and B with high values of hq and m lead, however, to relatively 
large values of the critical chemical potential /i qm at which the transition from nuclear to quark 
matter occurs. One also finds large values of /i e s,qm and the associated critical baryon density 
n qm in the quark matter phase. This sheds doubts on the reliability of this computation. One 
may argue that for such high values of /x e g there is no good reason why a separate minimum 
for nuclear matter should persist. The prediction of a first order transition for the saturation 
scenario (C) seems more robust in this respect. 

We next present a short description of the dominant effects that lead to our picture of a first 
order quark-hadron phase transition. This should also give an impression of the substantial 
uncertainties still inherent in this picture. The dominant mechanism for a possible first-order 
quark-hadron phase transition is the rapid decrease of at a — due to the fluctuations 
of massless quarks, whereas at the potential minimum which corresponds to nuclear matter 
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the effect of the quark fluctuations is reduced by their effective mass and by qn > 0. For 
\i increasing beyond /i nuc the density of nuclear matter increases beyond n nuc and a^ nuc \ii) 
decreases (scenarios A, B). In our crude picture this continues until fi e g reaches the value 
AV(<3 r ^ nuc ' ) ) (cf. eq. (51)). At the corresponding density the equation of state becomes very stiff, 
similar to the saturation scenario C discussed in section 6 (eq. (55)). The density can further 
increase because of quark contributions only once /i e s becomes larger than / u q (a (nuc - ) ). (For the 
parameter set C corresponding to the saturation scenario /j, e g must first exceed /^(ow) = 
334 Me V (/i > 383 Me V) before the density can increase beyond nuclear density.) As a result 
of the "frozen density" the rate of decrease of is also frozen according to eq. (3). On the 
other hand, the quarks always fully contribute to dU jd\i at a = (in the absence of current 
quark masses). For fj, ef[ > g H and ,u cff > /x n f(a nuc ) the potential at a = decreases therefore 
faster with fx than for the nuclear matter phase at a( nuc ) (cf. eqs.(17), (18) with q u (a = 0) = 
and g H (onuc) = qu). 

In the vicinity of a = the effect of the current quark masses should be included for a 
quantitative calculation. They push the minimum to positive a such that large enough quark 
masses typically destroy a possible first-order transition. In fact, at sufficiently high fi the 
effective potential (22) always has a new minimum near the origin at 



^ = 7^tf° ' (61) 



Here the mass parameter 
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CT=0 

m 2 = 2a 2 (A - 73 + 74 - 7s) -™l = + (74 - ^7s - ^75) cr 2 - 3ml . (63) 

corresponds to the curvature of the potential at the origin and m\o§ reflects the linear source 
term. At this minimum the effective quark mass m^" 1 -* = ha^" 1 ^ = haorn^ / m^/j) , vanishes in 
the chiral limit m 7r — * and for /i — > 00. For small enough m^ 111 - 1 we identify the corresponding 
phase with quark matter. For a vanishing current quark mass (j = 0, m 7T = 0) chiral symmetry 
is restored in this phase. In case of a first order phase transition and, in particular, for small 
current quark masses, one typically finds a situation where two different local minima at 
^(qm) an( j ^(nuc) coex j s ^ \ s ^ increases, the height of the potential for the quark matter 
phase U( qm )(fj,) = U(a^ qm \fi); /i) decreases faster than the one for the nuclear matter phase 
£y-(nuc)^ _ [/(^("^(/i); yu), where we remind that decreasing U corresponds to increasing the 
pressure p = —U. This can be seen directly from eqs. (19), (20), since cf( qm ) < cf( nuc ) and 
^(a;/ji) = 0. One concludes that for large enough \i the absolute minimum of U is always 
given by eqs. (61)-(63). 

Away from the chiral limit the quark-hadron phase transition is not characterized by a 
change of symmetry in our model 13 . It could therefore be of first order or a crossover. (A 
second order transition would require an additional tuning of parameters.) Our numerical 



13 As mentioned in the introduction, we do not take into account in the present approach the possible 
spontaneous breaking of color at high density [9] or in the vacuum [10], [11]. 
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evaluation of U for the three scenarios A, B, C shows that the existence of a first-order transition 
depends on assumptions about qu{a). For many reasonable functional forms we find indeed 
a first-order transition. In view of the remaining uncertainties it seems useful to establish 
general criteria for the occurrence of a first-order transition within our computation. Indeed, 
a first-order transition is guaranteed if a range of ji exists for which m^{ji) is positive and 
substantially larger than whereas a^ nuc \fj) remains of the same order of magnitude as vq. 
In this case one has tf( qn %) < a^ uc \/i) and the mass term at cr^ qm ^(/x) is well approximated 
by m^/i) and therefore positive. By definition the mass term at a^ nuc \ii) is also positive. 
Two local minima of U coexist for this range of /i. As /x is increased further the mass term 
at (f( qm )(/j) monotonically grows (cf. eq. (62)) thus excluding a crossover. Typical values for 
m from an extrapolation of the polynomial potential (28) for the parameter sets (A, B, C) 
are (833.4, 709.7, 584.9) MeV. For these values a first order transition would be guaranteed for 
m (pL eS ) k, 400 MeV or fe 9 26 ^ cV > (13.10,11.82,9.92) if <f (nuc) (/ieff) remains of order a . We 
use here a vacuum constituent quark mass of 330 MeV or a Yukawa coupling h q = 7.1. 

In summary, we infer a first-order transition if <5 r(nuc ' ) (/i e ff ) remains of order a for /x eff = 
(480, 433, 363) MeV for the models A, B, C. For low enough values of m (as, for instance, in 
scenario C) the value of /x e ff is low enough such that ^ nuc ^(/J e ff) is not expected to be much 
smaller than the value in nuclear matter at low pressure, as given in table 2. A first-order 
transition occurs then independently of other details of the potential. On the other hand, for 
large values of /j e fr the dependence of o r( - nuc - ) (/i e fr) is much more difficult to assess. It depends 
crucially on the way how the quarks are "switched on" , as expressed in the present formalism 
by the functional form of qn (&)■ A crossover or even a rather smooth transition become possible 
as well. 

Actually, a very natural scenario seems to be a first order transition at a critical value /x e ff, qm 
which is lower than /i q (cx ( ' imc * ) ) . In this case the quarks do not contribute in the nuclear matter 
phase and nucleons are absent in the quark matter phase. In the following we concentrate 
on this scenario which can be realized for the ansatz (26) with reasonable values of a c and 
(Tfj. Typical values of // e ff,qm f° r this situation are somewhat above ^h, say, // e ff,qm — (300 — 
400) MeV. The baryon density in the quark phase at this transition would be around three 
times nuclear density. These values occur naturally for values of /io somewhat below /x e fr,qm- An 
investigation of the coarse grained effective potential in the framework of the average action 
for a nonvanishing baryon chemical potential [40] finds values of /io only slightly above the 
constituent quark mass. This can be interpreted as an information about the potential Uq near 
the origin and supports the above scenario. 

In order to estimate the critical value // e ff,qm for the quark hadron phase transition in this 
scenario we equate the pressure in the quark matter phase (for the approximation m % = 0) 



Here and p are the nucleon mass and pressure corresponding to /j, e g = // n f, respectively. 
For qii not much larger than g nuc one may neglect p. Inserting two typical sets of values, 




(64) 



with the one in the nuclear matter phase 




(65) 
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Ho = 320 MeV, q H = 1.05(1. 2)g nuc , m N = 0.7(0.6)m N , one obtains /J eff)qm = 390(440) MeV. 
This corresponds to a critical baryon density in the quark matter phase 

n qm = 3.4(4.9)n nuc . (66) 

For the scenarios (A) and (B) the existence of two minima found in our computation 
may well be an artefact of our inaccurate treatment of the transition from quarks to nucleons. 
Indeed, there are reasonable forms of qui®) for which the nuclear matter minimum has reached 
small values of <7 nU c already for substantially smaller /j, e s . This would favor a smooth transition. 
One should remember, though, that the estimate of fics depends crucially on the value of m in 
eq. (62), which is only poorly known. In fact, the observed meson masses and decays contain 
only very limited information about the behavior of Uq near the origin. We do not expect a 
polynomial expansion of Uq around o"o to lead to a very good approximation of the potential in 
the vicinity of the origin. It is certainly possible to extrapolate a form of Uq which is compatible 
with nuclear physics constraints in the region 0.6<7o < a < 1.5<j to a wide range of parameters 
/io and m characterizing the behavior of Uq near o = 0. Furthermore, a possible cr-dependence 
of or M w would substantially affect the ratio Hcs/h- I 11 particular, a smaller value of 
g( u )/M u for the quarks (near a = 0) would enhance the effective chemical potential for given 
H in the quark phase, thereby shifting the transition to lower values of \x. We conclude that 
the spread in the values in table 3 (especially those corresponding to sets A and B) should be 
considered as an illustration of the uncertainties still inherent in the polynomial extrapolation 
rather than as actual predictions (which we expect closer to eq. (66) in case of a first order 
transition). This uncertainty is reduced significantly once independent information about the 
behavior of Uq near a = becomes available as, for example, from ref. [40]. 

We have also computed the surface tension £ qm for the quark-hadron transition at the 
critical /i qm . It turns out to be much larger than the one between the nucleon gas and nu- 
clear matter. The quantitative value is given in table 3. The surface tension depends, however, 
strongly on the details of the transition from quark to nucleon degrees of freedom (e.g., a c 
and oh). Stability of nuclear matter requires the critical chemical potential /i qm for a possible 
quark-hadron transition to be above /z nuc as realized for our parameters. At the quark hadron 
phase transition the quark mass M q qm in quark matter is much smaller than in nuclear matter. 
Nevertheless, it is substantially larger than the current quark mass. We quote the value of the 
order parameter cx qm = ^ qm ^(// qm ) for the quark phase in table 3 together with the correspond- 
ing quark mass. For // = /j qm the values in the nuclear matter phase are a^ nuc ^ (/x qm ) — 24 MeV, 
Mg ~ (165 - 170) MeV. Since // qm may exceed the effective strange quark mass, the strange 
quarks could play a role for this transition in real QCD. 

To summarize this section, our first computation exhibits a first-order phase transition 
between nuclear and quark matter at a critical density which is a few times nuclear density. If 
such a transition really occurs, our values for the density seem much more realistic than the 
very low values typically obtained in simple quark model computations. This underlines the 
importance of the correct treatment of the baryons in the nuclear matter phase. We have also 
seen, however, that the uncertainties remain very substantial and a smooth behavior remains 
also conceivable. The unknowns we have encountered will be present in any realistic mean 
field-type treatment of the transition between quark and nuclear matter. Whereas high density 
quark matter can perhaps be dealt with rather reliably at high enough baryon density, the main 
problem concerns the behavior of nuclear matter at high density for which confinement effects 
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(the binding to baryons) cannot be neglected. Any discussion of a phase transition needs 
knowledge about both phases concerned. The high density nuclear matter phase therefore 
needs always to be understood quantitatively and a pure quark model cannot give a reliable 
description (unless the "binding of baryons" is somehow incorporated). This raises substantial 
doubts about the applicability of many mean field statements about this transition in simple 
quark models. 

9 Conclusions 

We have presented here a new method for the computation of the dependence of the free 
energy on the chemical potential. It is based on an approximate solution to an exact functional 
differential equation. This method allows us to put mean field theory into a more systematic 
context. Chiral symmetry is explicitly implemented and phenomeno logical information about 
pion masses and decay constants is taken into account. Expressions which are close to mean field 
theory describe the difference in the free energy between vanishing and no n- vanishing chemical 
potential. They can be considered as the leading order in a series of systematic truncations of 
the exact differential equation (6). On the other hand, the free energy for vanishing chemical 
potential is not reliably described by mean field theory. Many relevant characteristics of this 
quantity can, however, be inferred from observation. 

Perhaps the most important new feature in our approach is that quark and nucleon fluctua- 
tions can be treated simultaneously within the same computation of the free energy. Only this 
allows a simultaneous description of the nuclear gas-liquid transition and the transition from 
nuclear to quark matter. A method which can deal both with nucleons and quarks is crucial 
for any quantitative treatment of a possible phase transition from nuclear to quark matter. 
Actually, for this transition the main difficulty lies in the understanding of the "low density 
phase". This phase is nuclear matter at a critical density of perhaps several times nuclear 
density, where standard nuclear physics is not of much help and the binding of quarks into 
nucleons nevertheless remains an important ingredient. Simple quark models not accounting 
for this binding, like NJL-type models, are insufficient for a description of this transition. 

We take here a very simple approximation where we use quark and nucleon degrees of 
freedom in their appropriate momentum ranges. For high momenta, q 2 > q\, the quark meson 
model gives a useful approximation. For small momenta, q 2 < q^, the effects of confinement 
have to be taken into account and we describe the carriers of baryon number as nucleons. Our 
simple model leads to a unified description of the nucleon gas, nuclear matter and transition 
to quark matter. The appearance of three phases of strongly interacting matter is related to 
three distinct minima of the effective potential for the a-field. (Typically only two coexist 
simultaneously.) This rich structure is a consequence of the fact that more energy levels fall 
below the Fermi energy for nucleons than for quarks. This results in a substantial enhancement 
of the density or, equivalently, the /^-dependence of the free energy due to low momentum 
nucleon fluctuations. The enhancement of the effect of a Fermi-gas of massive nucleons as 
compared to a gas of quarks with constituent masses only about a third of the nucleon mass 
is actually a huge factor of 27. This large enhancement relies only on the different effective 
masses and, therefore, Fermi surfaces of nucleons and quarks. Consequently, this property is 
rather independent of other more detailed features of the model. 
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The use of nucleoli degrees of freedom and the corresponding large density enhancement 
factor as compared to quark degrees of freedom is crucial for a any realistic description of the 
liquid-gas nuclear transition. It also shifts the transition from nuclear to quark matter to a 
larger chemical potential and baryon density, as compared to a description of the fermionic 
fluctuations in terms of quarks alone. This results in a substantial pressure at the coexistence 
between quark and nuclear matter. Standard Nambu-Jona-Lasinio type quark models fail to 
take into account this "density enhancement" and, therefore, typically predict a transition to 
quark matter at too low densities. 

The gas-liquid nuclear transition can be described within the present approach in a quan- 
titative way, at least for low temperature. We present here two classes of models, one similar 
to the a — w-model widely used in nuclear physics, the other a "saturation model" where the 
short distance repulsion between nucleons arises not only from vector-meson exchange but 
also reflects the transition to quark degrees of freedom at high density. For a suitable choice of 
parameters of the effective meson masses and interactions in vacuum our models give a success- 
ful description of the nuclear droplet model, consistent with the observed nuclear density, the 
binding energy per nucleon, the compression modulus and the nucleon mass in nuclear matter. 
They explain why nuclear density is approximately independent of the baryon number of a 
nucleus. We also obtain a realistic value for the nuclear surface tension and we have computed 
small corrections to the baryon density and the mass formula for nuclei due to nonvanishing 
pressure. Isospin violation and electromagnetism can be incorporated easily in our model. This 
should give a reliable equation of state for neutron stars in the region of moderate densities. 
Our model predicts coupling constants which directly enter the effective nucleon-nucleon po- 
tential. Comparison with nucleon scattering experiments will provide an interesting test in the 
future. 

A potential shortcoming of these models are the large meson self-interactions in vacuum 
which are needed for a realistic description. Within a linear a-model they correspond to large 
coefficients of terms ~ <r 6 , <r 8 or <t 10 which have not been found in previous renormalization 
group studies. It remains to be seen if a more complex structure of vacuum expectation values 
in the nuclear matter phase and in the vacuum - including spontaneous breaking of color - 
can ease this problem. The generalization of our approach to additional order parameters is 
straightforward. One may consider our quantitative results as a "prototype calculation" for 
understanding the properties of nuclei from QCD. It can be adapted to more complex settings. 
Particularly important will be a reliable computation of the vacuum-effective potential. 

Within our approximations we find a first order phase transition between nuclear and quark 
matter at high density and vanishing temperature. The critical density is typically three to 
five times nuclear matter density. We emphasize, however, that some important information is 
still missing for a quantitative understanding of the quark-hadron transition: The first problem 
concerns the appropriate formulation of a coarse grained effective potential and a determination 
of its shape for the vacuum. This is needed since first-order transitions require a non-convex 
coarse-grained potential whereas the inclusion of long wavelength fluctuations leads to a convex 
effective potential. Within the linear quark meson model we have addressed this issue in the 
context of the average action [40] . The second loose end is a more detailed understanding of the 
change from quark to baryon effective degrees of freedom. This concerns primarily the behavior 
of nuclear matter at densities much larger than nuclear density. It is therefore very relevant for 
a quantitative description of the quark-hadron transition. Furthermore, the neglected strange 
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quarks could play a relevant role at sufficiently high density. Finally, the presence of additional 
"color symmetry breaking" expectation values are a crucial ingredient for the understanding 
of states with high baryon density. Again, this could influence the nature of the transition 
from nuclear to quark matter [41]. In view of all these uncertainties we conclude that no 
definite statements about the nature of this transition can be made at present. Both a genuine 
first-order phase transition and a relatively smooth change remain possible. The uncertainties 
mentioned here are common for other analytical approaches and constitute a major difficulty 
for a quantitative understanding of the quark-hadron phase transition. 
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